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Abstract. Although several models of f(R) theories of gravity within the Palatini 
approach have been studied already, the interest was concentrated on those that have 
an effect on the late-time evolution of the universe, by the inclusion for example 
of terms inversely proportional to the scalar curvature in the gravitational action. 
However, additional positive powers of the curvature also provide interesting early-time 
phenomenology, like inflation, and the presence of such terms in the action is equally, 
if not more, probable. In the present paper models with both additional positive and 
negative powers of the scalar curvature are studied. Their effect on the evolution of 
the universe is investigated for all cosmological eras, and various constraints are put 
on the extra terms in the actions. Additionally, we examine the extent to which the 
new terms in positive powers affect the late-time evolution of the universe and the 
related observables, which also determines our ability to probe their presence in the 
gravitational action. 


f sotiriou@sissa.it 


Constraining f(R) gravity in the Palatini formalism 2 

1. Introduction 

General Relativity has been an extremely successful theory of gravity for almost 100 
years now. It has passed all tests on scales relevant for the solar system. However, all 
of these tests are actually relevant to the Post-Newtonian regime and not to the full 
version of the theory. At the same time, during the last decades, several indications 
have appeared, leading to the thought that maybe General Relativity is the effective 
part of a more general theory of gravitation, applicable only on solar system scales. 
Some of these indications are of theoretical origin and as an example one could name 
the difficulties in finding the quantum counterpart of general relativity and the fact 
that it has not yet been possible to include it in any grand unification scheme. The 
others are of observational origin, like the fact that relativity cannot explain the flat 
rotation curves of galaxies without introducing unseen dark matter and, of course, the 
fact that it cannot explain the current accelerated expansion of the universe without 
the introduction of both dark matter and dark energy. Also, in order to explain these 
cosmological observations, one has not only to introduce these exotic componets, but 
also to assume that they sum up to 96% percent of the total energy content of the 
universe. This, of course, does not sound appealing at all, since the nature of these 
components still remains a mystery, and this has lead to several attempts to modify or 
extend General Relativity in different ways, addressing the dark matter [3 El El El (3 EJ 
or the dark energy 0 El 0 EH El problem. 

An alternative route is to modify General Relativity by abandoning the simplicity 
assumption that the action should be linear in the scalar curvature R. So, instead of 
using the Einstein-Hilbert Lagrangian one can use a more general one, which depends 
on a generic function f(R) (see (I2( HB1 IHj for early studies). Hence, such theories 
are called f(R) theories of gravity. Such modifications are not completely arbitrary as 
higher order terms in curvature invariants seem to be present in the effective Lagrangian 
of the gravitational field when quantum corrections or string/M-theory are considered 
Ha naira nans EDi. Therefore, /(/?,) theories seem also to address the theoretical 
concerns about General Relativity that were mentioned earlier. Phenomenological 
interest up until recently was restricted to the inflationary behaviour exhibited when 
positive powers of R are present in the Lagrangian EH- However, it was seriously 
increased when it was shown that a term inversely proportional to R can lead to late 
time expansion 1221 EBl EU- Such theories, appealing as they may be, are unfortunately 
not free of problems. First of all they lead to fourth order differential equations 
which are difficult to attack. Additionally, it is doubtful whether they can pass the 
known solar system tests [23 EE! or whether they have the correct Newtonian limit 
H3EIIE3EH1- It is possible that more sophisticated models may exhibit behaviour 
closer to the expected one ra, but this requires significant fine tuning of the various 
parameters. The most important problem of these models, however, is that they lead 
to unavoidable instabilities within matter in a weak gravity regime jP> 1 j . 

There is a further modification of gravity that one can consider which does not 
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necessarily involve any modification of the action, but rather the use of a different 
variational principle. This variational principle, known as the Palatini formalism 
due to a historical misconception, since it was not introduced by Palatini but by 
Einstein himself, treats the metric and the affine connection as independent geometrical 
quantities. One has to vary the Lagrangian with respect to both of them to derive the 
field equations, in contrast with the standard metric variation, where the Lagrangian 
is varied with respect to the metric alone, and the connections are assumed to be the 
Christoffel symbols of this metric. It is also known that when the Einstein-Hilbert action 
is used, the Palatini variational principle leads to the Einstein equations, just like the 
standard metric variation This is not true, however, for a more general action. 

When used together with an f(R) Lagrangian, the Palatini formalism leads to 
second order differential equations instead of the fourth order ones that one gets with 
the metric variation. At the same time, in vacuum, they straightforwardly reduce to 
standard General Relativity plus a cosmological constant J53]. This ensures us that, 
firstly, the theory passes the solar system tests, and secondly, that interesting aspects 
of GR like static black holes and gravitational waves are still present. It is also free of 
the instabilities discovered in inn for f(R) gravity in the metric approach, which arise 
within matter in a weak gravity regime. Even though there was initially some debate 
concerning the Newtonian limit 551 . a recent paper |2S] seems to have settled this 
matter, showing that these theories indeed have the correct behaviour. One more reason 
to study such a metric-affine variation is theoretical completeness. Both approaches give 
exactly the same result when applied to the Einstein-Hilbert action. Thus there is no 
real criterion so far about which one of them is better to use. Additionally, though, 
the Palatini variation seems to be more general since it yields GR without the need to 
pre-specify the relation between the metric and the connections. Finally it should be 
mentioned that the Palatini formalism is also closer to the picture of a Hilbert space, 
since is assumes that there are two sets of independent variables, and therefore seems 
more appealing when trying to quantize gravity. 

In jHE], Vollick showed that using the action of 12a, which includes an R 1 term, 
together with the Palatini variational principle can lead to a theory of gravity that 
predicts late time accelerated expansion of the universe. Several models with similar 
behaviour followed (for a short review see IHZj and references therein). Later it was 
proved that including an R 2 term in the action could not drive inflation J3H within this 
formalism. However, in 22j it was shown that this is not true for higher order terms 
in general, but merely that R 2 constitutes an exception due to the construction of the 
theory, and a model was presented in which both an R 3 term and an R~ l term were 
present in the action. It was demonstrated that such a model can account both for early 
time inflation and late time accelerated expansion. 

It is important, of course, to go beyond the qualitative results and use the numerous 
observations uniimmiii to get quantitative ones. Such a study was performed in 
jPHh Assuming that the gravitational action includes, besides the standard linear term, 
a term inversely proportional to R , the authors used four different sets of cosmological 
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data to constrain it. These are the Supernovae type la gold set [U3, the CMBR shift 
parameter [SS], the baryon oscillation length scale [4B] and the linear growth factor 
at the 2dF Galaxy Redshift Survey effective redshift ga eb]. However, as stated in 
the conclusions of the restricted form of f(R), including only a term inversely 
proportional to i?, prevents the study from being exhaustive. 

Even though several models of f(R) gravity in the Palatini formalism have been 
studied, most interest was concentrated on those having terms inversely proportional to 
the scalar curvature. Bearing in mind that including an additional term proportional to 
a positive power of the curvature in the action leads to interesting phenomenology |39| 
and of course is equally, if not more, reasonable to have present in the action, here the 
behaviour of models with both positive and negative powers of R will be investigated 
(we will be using the term positive powers of R to imply positive powers additional to 
R itself). By generalizing the results of [22] it will be shown that such models lead to 
curvature driven inflation and late time cosmic expansion. However, a more detailed 
study reveals that such a curvature driven inflation is not interesting but is actually 
problematic, because there is no way to put an end to it, and should therefore be 
avoided. Another important question is how to constrain terms in positive powers of 
R. This problem will be addressed here in detail. One of the possible constraints will 
arise if we want to avoid the undesirable curvature driven inflation mentioned before, 
but also more standard ones, related to the later evolution and the Newtonian limit 
will be considered. The results will also be used to check whether the presence of such 
terms could affect the results of M- Additionally, a thorough discussion will take place, 
which can be considered relevant, not only to the results of this paper, but also to studies 
similar to H- 

The rest of the paper is organized as follows: In section [21 the Palatini formalism is 
briefly reviewed. In section El a model including both negative and positive powers of the 
scalar curvature in the gravitational action is examined. Its behaviour is investigated 
for different cosmological eras and the possible constraints for the positive power term 
coming from different origins are obtained. Section |U contains a thorough discussion of 
the physical consequences of the derived results. Also, the results of J0] are interpreted 
and discussed and we examine, using the results of section 0 whether they can actually 
be consider general enough to include gravitational Lagrangians with positive powers in 
the scalar curvature as well. Section El contains conclusions. 

2. The Palatini formalism 

We start by very briefly reviewing the Palatini formalism for a generalized action of the 
form 

s = hJ d ' x ^^ R ) +^m, (i) 

where k = 8ttG and Sm is the matter action. For a detailed study see [3S]. 
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As mentioned in the introduction, within the Palatini formalism, one treats the 
metric and the affine connections as independent quantities. Therefore, we have to vary 
the action with respect to both of them in order to get the field equations. Varying with 
respect to the metric g liu gives 

f\R)R, v - l -f{R)g„ = k (2) 


where the prime denotes differentiation with respect to R and the stress-energy tensor 
Tfu, is given by 


, _ 2 SS M 

^ V=9Sg^ 


(3) 


Varying the action with respect to T^ IU and contracting we get 


Va\f(R)y/=g<r] = 0 , 


(4) 


where V denotes the covariant derivative with respect to the connections. It is 
straightforward now that the connections are the Christoffel symbols of the conformal 
metric h^ u = f'{R)g tu , (see [3BJ). Thus, we get 

+ fj-, [Kw - . (5) 

where {^} denotes the Christoffel symbol of g tw . Finally, by contracting eq. (0) one 
gets 


f{R)R-2f{R)=nT. (6) 

Eq. (0) can in general be solved algebraically for R = R{T ), if of course the functional 
form of / allows it. 

If we are interested in cosmological solutions we can consider the spatially flat 
FLRW metric, 

ds 2 = — dt 2 + a[t) 2 8ijdx l dxP (7) 


and the perfect fluid energy-momentum tensor Tf 

Friedmann equation is given by 
/ . \ 2 


H + 


1 r 

2 r 


1 k(p + 3 p) lf_ 
6 f 6 f 


diag(— p, p,p, p). The generalized 


( 8 ) 


We can use eq. (0 together with the conservation of energy to express R as a function 
of R (see also Iffifl): 


3 H(Rf - 2/) 

Rf" ~ f 


(9) 


Using eq. © to re-express /'(= f"R) and assuming that the universe is filled with dust 
(p = 0) and radiation (p = p/3), after some mathematical manipulation, eq. © gives 

1 2 np + Rf'-f 


H 2 = 


6/' (I _ 3 /"(fl/'-2/) \ 

w 2 ) 


2 ’ 


and p = p m + Pri where p m denotes the energy density of dust and p r the energy density 
of radiation. 
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What we would like to investigate here is the cosmological behaviour of a model with 
both positive and negative powers of the scalar curvature in the gravitational action. 
Such a general study seems to be a tedious analytical task given the complexity of the 
functions involved and the non-linearity of the equations. An alternative could be to do 
a numerical analysis which would however lead to lack of generality, and therefore should 
not be so appealing. Fortunately, there seems to be a more elegant way to approach this 
problem, following the lines of [32]. In some part of this paper we will leave the function 
/ unspecified and try to derive results independent of its form. However, in some cases 
we will also adopt the following representation for /, suitable for our purposes: 


f(R) = ^R d + R 


6+1 
£2_ 

R b 


( 11 ) 


with ci, e 2 > 0, d > 1 and b > 0; b = 0 corresponds to the ACDM model when e\ —> oo. 
The dimensions of e\ and e 2 are (eV) 2 . Part of our task in this section will be to 
constrain the values of e±. For the value of e 2 however, no extended discussion is really 
necessary, since we already know that in order for a model to be able to lead to late time 
accelerated expansion consistent with the current observations e 2 should by roughly of 
the order of 10~ 67 (eV) 2 . 

In [22] the specific case with 6 = 1 and d = 3 was discussed and it was shown 
that such a model can account for both early time inflation and late time accelerated 
expansion, and also give a smooth passage from one to the other during which standard 
cosmological eras can occur. Here we generalize this approach. First of all, let us see 
how a model with a general function / would behave in vacuum, or in any other case 
that T = 0 (radiation, etc.). If we define 

F(R) = f\R)R — 2/(i?), (12) 

then eq. © reads 


F(R) = kT, 
and for T = 0, 

F(R) = 0. 


(13) 

(14) 


Eq. (11411 is an algebraic equation which, in general, will have a number of roots, R n . Our 
notation implies that R is positive so here we will consider the positive solutions (i.e. 
the positive roots). Each of these solutions corresponds to a de Sitter expansion, since 
R is constant. If one wants to explain the late time accelerated expansion one of these 
solutions, say i? 2 , will have to be small. If in addition to this we also want our model 
to drive an early time inflation, there should be a second solution, Ri corresponding to 
a larger value of R. 

For example, introducing in eq. © the ansatz given for / in eq. dm), one gets 

^z^R d+b - R b+1 + (6 + 2)c2 +1 = 0. (15) 

e i 
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Ri ~ Ci, i ?2 ~ ^ 2 - (16) 

These solutions can act as seeds for a de Sitter expansion, since the expansion rate of 
the de Sitter universe scales like the square root of the scalar curvature. Notice that 
the rest of the solutions of eq. (fUj) , for R > Ri and R < Ro will not be relevant here. 
During the evolution we don’t expect, as it will become even more obvious later on, 
that R will exceed i ? 2 or become smaller than R±. 


3.1. Early times 


Following the lines of we can make the following observation. Since we expect 
the matter to be fully relativistic at very early times, T = 0 and consequently R is 
constant. This implies that the second term on the left hand side of eq. (JBJ) vanishes. 
Additionally, conservation of energy requires that the first term on the right hand side of 
the same equation scales like a{t )~ 4 , which means that the second term on the same side, 
depending only on the constant curvature, will soon dominate if f(R) is large enough 
for this to happen before matter becomes non relativistic. R can either be equal to R\ 
or to i? 2 - Since we want Ro to be the value that will provide the late time acceleration, / 
should be chosen in such a way that /(/?, 2 ) will become dominant only at late times when 
the energy densities of both matter and radiation have dropped significantly. Therefore, 
if we want to have an early inflationary era we have to choose the larger solution R \, 
and / should have a form that allows f(Ri) to dominate with respect to radiation at 
very early times. The Hubble parameter will then be given by 


H ~ 


Gf'(RiV 


(17) 


As an example we can use the ansatz given in eq. (ED- The modified Friedmann equation 
reads 


H ~ \Js(d+ 1 )’ (18) 

and the universe undergoes a de Sitter expansion which can account for the early time 
inflation. 

Sooner or later this inflationary expansion will lead to a decrease of the temperature 
and some portion of the matter will become non relativistic. This straightforwardly 
implies that R will stop being constant and will have to evolve. f'(R ) plays the role 
of the conformal factor relating the two metrics g /lu and h llu , and therefore we do not 
consider sign change to be feasible throughout the evolution of the universe. We also 
know that in a certain range of values of R it should be close to one. This is the case 
because there should be a range of values of A, for which f(R) behaves essentially like 
R , i.e. our theory should drop to standard General Relativity, in order for us to be able 
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to derive the correct Newtonian limit (see H and section Id.dll . Together with T < 0, 
the above implies the following: 

f\R) >0, F < 0, V R 2 < R < Ri- (19) 


Since F is a continuous function keeping the same sign in this interval and F(Rf) = 
F(R 2 ) = 0, there should be a value for R, say R e , where F'(R e )= 0, i.e. an extremum. 
Eq. fTilj) implies that the time evolution of R is given by 


R = kT/F\R). 


( 20 ) 


Differentiating eq. m we get 


F’ = f"R-f. (21) 

Using the fact that f = f'R, and using eq. m to express f" in terms of F', f and R 
one can easily show that 


t 

r 


F'+r 

Rf 


R. 


( 22 ) 


The constraints given in eq. m imply that for R e < R < R\, F' > 0. An easy way to 
understand this is to remember that F is negative in that interval but zero at R = R\ 
and so it should be an increasing function (see also fig. HJ). Since, f and R are also 
positive, then what determines the sign of /'/ f in the neighbourhood of R e is the sign 
of R. 



Figure 1 . The behaviour of a general function F(R) over the interval i?i > R > R 2 . 
R e denotes the value of R where F has a minimum. From this graph one can easily 
see that F' is possitive when R e < R < R\ and negative when R 2 < R < Re- 


Let us see what will happen if we require R to decrease, i.e. R < 0. Eq. © implies 
that as R —> R e , R —> —oo if T ^ 0, since F'(R e ) = 0. Therefore, f'/f —> —oo and 
using eq. © we can infer that F[ — > oo. Physically, the above imply the following: R 
has no way to decrease to a value less than R e without giving the universe an infinite 
expansion. In practice, any attempt for R to approach R e would lead to a dramatically 
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fast expansion, np until non-relativistic matter fully dilutes, and R settles back to 
Ri. Thus, once the curvature terms in the modified Friedmann equation dominate the 
evolution, there is no turning back to matter domination through a continuous process. 
The two vacuum solutions Ri and R 2 seem to be somehow disconnected in the evolution, 
and R has to remain in the region close to only one of them. This is a general statement 
independent of the form of matter that is present, since T was left unspecified in its 
derivation. So, even though, as shown in j22], including positive powers of R in the 
action can lead to early time inflation, there seems to be no graceful exit from it. The 
only alternative left would be to consider that due to some other physical and non- 
classical process, the validity of the equation presented here ceases to hold for some 
time interval, which, however seems highly unphysical. 

Since it seems impossible to provide an exit from this gravity driven inflation it 
seems reasonable to check if we can at least totally avoid it. If we choose as our initial 
solution R 2 instead of Ri then the curvature terms will not dominate as long as R is 
constant, ffowever, there is still one subtle point. At some stage during the evolution 
the energy density of non-relativistic matter will have to rise sooner or later, forcing 
R to change its value. It is also reasonable to assume that if inflation is not driven 
by curvature we will have to adopt a more standard approach to guarantee that it will 
happen, like an inflaton field. It is obvious, however, keeping in mind the previous 
discussion, that one would want R to always be less than R e and this will impose a 
constraint which will depend on the functional form of /. For example if one assumes 
that / is described by the ansatz given in eq. ED, then, considering ordinary matter, 
R < R e at all times implies that ^ np m at all times. Let us also consider the case 
of a slow-rolling inflaton field, <f>. Then, if we denote its energy density by p^, and its 
pressure by p^, we have, as usual, 


P0 = ^0 2 + ^(0), 

(23) 

p<p = - V((/>), 

(24) 


where 1/(0) is the scalar field potential. During the period when 0 dominates the 
evolution T = 0 2 — 41/(0), and since slow-roll implies that 0 2 <C 1/(0), T « —41/(0). 
Therefore, if we want R < R ei so that inflation proceeds as usual, then 1/(0) at 
all times. 

If R is less than R e for all values of T then it is easy to verify that everything 
will evolve naturally after the end of inflation. For R 2 < R < R e , F' < 0 and as non- 
relativistic matter dilutes, p m < 0, so from eq. (ED we see that R < 0 so that R will 
have to decrease to reach the value R 2 asymptotically. 

The above discussion, is not relevant of course if the only positive power present 
in the action is R 2 . In this case, due to the form of F , this term does not appear in 
eq. m (see also [2SJ). This specific case has been studied in E2J. One thing that is 
worth commenting upon, before closing this discussion, is the following. The constraint 
f > 0 (see eq. (PSD), which is implied by the fact that f plays the role of the conformal 
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factor relating the metrics and h tw , can, depending on the form of / impose a further 

constraint on the value of the constants in front of the positive power terms. In 
inflation driven by an inflaton held was studied in the presence of an R 2 term. The 
authors derived a constraint for the constant appearing in front of the R 2 term in the 
action, by requiring that the square of the Hubble parameter is positive during the 
kinetic dominated phase. This constraint is exactly what would one derive by requiring 
f to be always positive. 


3.2. Big bang nucleosynthesis 


Let us now turn our attention to the following cosmological era, radiation domination 
and big bang nucleosynthesis (BBN). Current observations indicate that the standard 
cosmological model can fit the data related to the primordial abundances of light 
elements. On the other hand, how a modified gravity model like the one discussed 
here would fit those data has not been yet worked out. However, there is little room for 
modifying the behaviour of the Friedmann equation during BBN and it seems reasonable 
to ask that the model under investigation should resemble standard cosmology during 
these standard cosmological eras Hi. This implies that eq. m should be similar to 
the standard Friedmann equation 

H 2 = -ftp. (25) 

By comparing eqs. m and ( 125 )) one see that during BBN 


1 

3 /W-2/) 
2 f'(Rf" ~ /') 
Rf - 3/ ~ 0 


(26) 

(27) 

(28) 


To make the picture clearer we give the explicit expressions for f and f" when / is 
given by eq. (ED: 


/' 


yd— 1 


,R~ " e 

d——— T 1 T b 


6+1 

2 


A -1 


Rb+V 


(29) 


p?d —2 ,6+1 

f" = d ( d ~ l)-^pr - b(b + l)-^2 ■ (30) 

Let us for the moment assume that the term inversely proportional to R is not present. 
In order for condition (El to be fulfilled e\ Rbbn ■ This is the natural constraint on 
the value of e± imposed when one asks for the model to have almost identical behaviour 
to the standard one during BBN. Once t\ is chosen to have a large enough value all three 
constraints El, El and EBP are easily fulfilled and the modified Friedmann equation 
(Cl becomes identical to the standard one, eq. El, for the relevant values of R. The 
above constraint can be viewed as a sufficient constraint for the model to be viable but 
not as a necessary one. However, one could also claim that, even if the modifications 
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in the Friedmann equation do not necessarily have to be negligible, they should at least 
lead to second order corrections and not affect the leading order. This implies that e\ 
should definitely be larger than Rbbn- 

We have, however, neglected the presence of the term inversely proportional to R. 
In the absence of positive powers this term should be negligible during BBN since Rbbn 
is much larger than e 2 . This picture may change if we consider the full version of the 
model. Eq. © can take the following form 

Rf-2f = -K P 0 m (l + z) 3 , (31) 

where p° m is the present value of the energy density of non-relativistic matter and z is 
the redshift. We have assumed here that a 0 = 1. Using eq. m one can derive how 
R will scale with the reshift. If e\ Rbbn, then for all the evolution of the universe 
after BBN, R scales almost like (1 + z) 3 . This indicates that, since BBN takes place 
at a very high redshift, Rbbn is indeed much larger than 62 - If however, one assumes 
that e± is large enough to alter the behaviour of eq. (EH), then R will have a milder 
scaling with the redshift, meaning that Rbbn can get very close to e 2 . Then the three 
constraints EH}, EH and EH might not be fulhlled not only due to the presence of the 
positive power of R , but also because of the presence of the negative one. This will be 
a secondary effect related to the term with a positive power as shown earlier. It will be 
avoided if again e x Rbbn and will be subdominant if e\ is just smaller than Rbbn- 
Unfortunately, since the value of Rbbn is very much model dependent it is difficult to 
turn this constraint into a numerical one. 

3.3. The Newtonian limit 

Before, examining the behaviour of the models under discussion at late times we will 
Erst check their Newtonian limit. The formalism for doing this has been derived in 
«2§]. Since we are using a metric-affine variational principle, the connections are not the 
Christoffel symbols of the metric g^, but of the conformal metric 

v = f(R)g^- (32) 

Expressed in terms of h^ the held equations read 

V'l = (33) 

It is easy to see that if f — 1, eq. EH reduces to the Einstein equation and the metrics 

and g /w coincide. Thus the model will have the correct Newtonian limit if f is very 
close to 1 for curvatures of the order of the solar system one. Current tests have an 
accuracy of 1CU 16 which implies that 

R ss < (34) 

where R ss is the scalar curvature for the solar system density. We can use eq.© 
to express the curvature in terms of the matter density. For the solar system p ~ 
10 - 11 g/cm 3 and after some manipulations we get 
€\ > l 0 16 /U - 1 U 45 (eV) 2 . 


Rfiu — -Rh^ + (f — 1 ) 


/? -A 

2 jo 


( 35 ) 
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In all the above we have neglected the contribution of the negative power of R since its 
contribution to f for this range of densities is of the order of 10 -21 ( b+1 b See also J22] 
for a more detailed discussion of the Newtonian limit. We also know that the current 
value of R , R 0 , is of the order of e 2 , so it easy to express eq. (E3) in terms of i? 0 : 

d > 10 16/(d ” 1)+21 i? 0 . (36) 

This condition will prove useful later. 


3 . 4 . Late times 


Now let us check what the behaviour of the modified Friedmann equation will be at 
late times. The scalar curvature R decreases with time to reach a value close to e 2 . 
Therefore the conditions (HID, (EZD and (|2%ll will at some point cease to hold because 
of the term involving the negative power of R. Any contribution of the term involving 
the positive power of R will be negligible for two reasons. Firstly, since the value of e\ 
should be such that these terms are already negligible during BBN, it is safe to assume 
that they will remain so all through the rest of the evolution of the universe. The 
same results can be inferred by using the constraints derived in section 13.11 Secondly 
there is also the constraint on the value of t\ coming from the Newtonian limit. As we 
mentioned earlier, the constraints coming from BBN are sufficient for a viable model 
but are not absolutely necessary, since it is not yet clear how well a model like that 
would fit the data from light element abundances. One could still assume that there 
is some slight contribution of the terms being discussed in the modified Friedmann 
equation during BBN, which in any case becomes even weaker at later times. The 
constraints coming from the early time behaviour are necessary but it is difficult to turn 
them into numerical ones. At the same time one can always claim that the early time 
evolution of the universe is not very well established and there might still be room for 
new physics there, affecting these constraints. However, the Newtonian constraint (USD 
is very straightforward and unquestionably necessary. This constraint will prove more 
than sufficient for our purposes in this section. 

The range of values of R which is of interest for late time observations is between 
the value of R at decoupling R dec and R 0 . Using eq. m it is easy to show that 
Rdec (1 + Zdec)~ 3 ~ Ro, or using as the value of the redshift at decoupling Zd ec = 1088 
one gets R dec ~ 10 10 i? 0 - Eq. (IHffll then implies that e x > 10 16 ^ d ~ 1 ' )+u R dec . Therefore 
the positive power term will be at least eleven orders of magnitude smaller than the 
linear term for any d > 2 and will become even more negligible as time passes. Thus, 
effectively 
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with extremely high accuracy for all times after decoupling. It is easy to see that the 
modified Friedmann equation of the model described in m will be identical to that of 
a model with no positive powers of the curvature {e\ —> oo) for late times. 

In |<3J the authors consider / to be of the form 


f(R) = R 1 + a 


R_ 

Hi 


/3-l N 


(40) 


where a and (3 are dimensionless parameters, with (3 < 1 (note that in our notation 
R is positive). This representation of the function / is very useful when one wants 
to constrain some dimensionless parameter. Comparing it with our ansatz, eq. m , 
we get d = S + 1, b = —f3 and e\ —> oo, since in eq. (SOI) there is no positive power 
of R. In order to constrain the values of a and (3 they use a rather extensive list of 
cosmological observations. The first quantity which they consider is the CMBR shift 
parameter [4511521 HH1 which in a spatially flat universe is given by 

rZdec dz 


n = 




H{z) 


(41) 


where Zd ec is the redshift at decoupling and Q m = np Q m / (37/q). When expressed in terms 
of the scalar curvature, eq. m becomes 

rZdec dz 


iz = \ a 
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f dR 


Jro {Rf 1 - 2/) 2/3 H(R)' 

Using the values for Zd ec and 1Z obtained with WMAP |1H] . namely Zd ec = 1088^2 and 
TZ = 1.716 ±0.062, they find that the best fit model is {a, (3) = (—8.4, —0.27). They also 
use the “Gold data set” of Supernovae BDj. What is important for this analysis is the 
expression of the luminosity distance, which in terms of R is 

dz 


D l {z) = (1 ± z) 
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Jro {Rf ~ 2/) 2/3 H{R )' 

Marginalizing over the Hubble parameter h, the authors again constrain a and (3 and 
the best fit model is {a, (3) = (-10.0,-0.51). Another independent observation which 
they use is that of the imprint of the primordial baryon-photon acoustic oscillations on 
the matter power spectrum. The dimensionless quantity A 


A = V^E{ Z1 ) 
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where E(z) = H(z)/H 0 can act as a “standard ruler”. The data from the Sloan Digital 
Sky Survey fTi provide a value for A, namely 


A = D v (z = 0.35) = 0.469 ± 0.017, 


0.35c 


where 


D v {z) = 


D m (z) 2 


cz 

W) 


1/3 


(45) 


(46) 


and Dm{z ) is the comoving angular diameter distance. The best fit model using this 
value is (a, (3) = (—1.1,0.57). Finally, in these three sets of data are combined to 
give a best fit for (a,/3) = (—3.6,0.09). 

The above observations are potentially very useful, of course, in studying the 
viability of a model like (IH . I will comment on the interesting results of (43, in the 
next section. However, as shown here, the modified Friedmann equation of a more 
general model like ED, which also includes positive powers of R , is effectively identical 
to that of d3DD at late times. Therefore, it is expected that the results of mi will remain 
unaffected by the inclusion of positive powers of the scalar curvature, since these terms 
have to satisfy the constraints derived in this section. 

Before closing this section 1 would like to briefly discuss the fourth scheme used in 
[33, in order to obtain constraints: large scale structure and growth of perturbations. 
The authors use the Jebsen-Birkoff theorem which, as they say, has not proved to 
hold for f(R) theories in general. Of course one would probably expect that it does 
hold, and once one accepts their assumption, their results come naturally. However, 
keeping in mind that these results cannot actually improve the constraints obtained 
with the three schemes already mentioned, there is no real reason to pursue such a 
study further. Therefore in the present paper 1 will not comment further or try to 
generalize these results. As the authors of mi correctly state, a more detailed analysis 
should be performed along the lines of E3- 


4. Discussion and physical interpretation of the constraints 

ft was shown that a model with a positive and a negative power of the scalar curvature 
added to the Einstein-Hilbert action can lead to the following evolutionary eras: a 
curvature driven inflation; a radiation dominated era in which BBN can take place, 
and during which the modified Friedmann equation can effectively be identical to the 
standard one, followed normally by a matter dominated era; a late-time era in which 
the universe undergoes a phase of accelerated expansion. This is in accordance with the 
results of j3J2j and requires certain constraints. 

However, as shown in section o a curvature driven inflation of this sort conies 
with serious problems. The curvature cannot drop to lower values in a continuous way 
without leading to an infinite value for the Hubble parameter. This implies that if the 
curvature is indeed dominant at early times, it is not possible to return to standard 
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cosmology at later times. Therefore, if higher order terms are indeed present in the 
action, one has to assume that the initial value of R is one that does not lead to early 
time curvature domination. Even in this case however, the theory comes together with 
a constraint. For a model described by eq. m this is €i p m , which should hold 
through all of the evolution of the universe. If, additionally, one asks for an inflaton 
field inflation to occur, €\ V (</>) at all times. These bounds imply that e j has to 
have a very large value, and are in accordance with the value of €\ expected from the 
perspective of effective field theory, ~ M 2 . Even though the analysis that leads to 
these result does not apply for the exceptional case of d — 2, as mentioned in section 
Id . 1 1 effectively the same constraints still have to hold, as shown in H- As mentioned 
here, they can easily be inferred from the positivity of f during the kinetic domination 
phase of the inflaton. The above results indicate that the presence of positive powers of 
R in the action cannot lead to interesting phenomenology and at the same time be in 
accordance with the established late time evolution of the universe. Such terms might 
be present in the action resulting from a more fundamental theory, but they would have 
to be very subdominant, even at early times and large curvatures. 

For the Friedmann equation not to be seriously modified during BBN we found 
that e 1 > Rbbn, which is a model dependent constraint. Finally requiring our model to 
have the correct Newtonian limit gave 6 \ > 10 16 /( d_1 )~ 45 (eV) 2 or in terms of the current 
value of the curvature on cosmological scales, e\ > 10 16// ( d-1 ) +21 -R 0 - This last constraint 
is of course a necessary one, whereas the first one is related to the phenomenological 
behaviour that our model should have. The constraint coming from BBN is definitely 
sufficient for the model to be viable and seems to be a necessary one, but a more detailed 
study is needed to reach a final conclusion about this. 

Making use of the derived constraints, and especially the one coming from the 
Newtonian limit, it has been shown that the presence of positive powers of the scalar 
curvature in the gravitational action cannot lead to any modification in the late time 
form of the Friedmann equation. Since this equation describes the evolution of the 
universe, the above statement implies that the late evolution of the universe is not 
affected by the positive powers of the scalar curvature present in the action. This can 
be rephrased in two interesting ways: The results of observational tests relevant to the 
late time evolution of the universe are insensitive to the inclusion of additional positive 
powers of R or observational tests relevant to the late time evolution cannot constrain 
the presence of additional positive powers of R in the gravitational action. The first 
expression makes it clear that the conclusions drawn from [41], or from any similar 
study, are actually more general than the authors claim. The second one however might 
be even more interesting since it implies that such tests are not sufficient to judge the 
overall form of the gravitational action. We will come to this shortly. 

We now discuss the results of [23j, which were reviewed in section El The best 
fit model for the combination of the different data sets suggests that their exponent /3 
is equal 0.09 (see eq. (HDD) and therefore favours the ACDM model, being well within 
the la contour. However, one gets different values for (3 when the different data sets 
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are considered individually. For the SNe data the best fit model has (3 = —0.51 and 
the baryon oscillations (3 = 0.57, both disfavouring the ACDM model, but also being 
mutually contradictory. The CMBR shift parameter gives (3 = —0.27 which again is 
significantly different from the other two values. Of course one might expect that the 
combination of the data will give the most trustworthy result. However, this is not 
necessarily true, and one could regard the discrepancies in the value of (3 coming from 
different observation as an indication that more accurate data are needed to derive 
any safe conclusion. It would also be interesting to study to what extent a model 
with (3 — — 1, which is the model commonly used in the literature, can individually 
fit the current data, and to compare the results with other models built to explain the 
current accelerated expansion, such as quintessence, scalar-tensor theories, etc. Bear 
in mind that the ACDM model has always been the best fit so far. However, the 
motivation for creating alternative models does not come from observations but from 
our inability to solve the theoretical problems that come up if one adopts the standard 
picture (coincidence problem, etc.). 

Notice that the results derived in section |21 hold also for 6 = 0 which is a model with 
a cosmological constant. This means that one could add to the statement of the authors: 
“data indicate that currently there is no compelling evidence for non-standard gravity”, 
that the data also indicate that currently there is no compelling evidence against non¬ 
standard gravity, since a model with 6 = 0 but d ^ 0 would give the exact same results as 
the ACDM model. This highlights the point made earlier, that the presence of positive 
powers of R in the action cannot be constrained using observational data relevant to 
the late time evolution. Such terms, however, could be seriously constrained using other 
types of data, such as the evolution of perturbations during inflation. 

5. Conclusions 

Starting by setting up the formalism needed, a model with both positive and negative 
powers of the scalar curvature, additional to the standard Einstein-Hilbert term in the 
gravitational action, was studied. The exponents and the coefficients of the extra terms 
were left arbitrary and all eras of cosmological evolution were studied. It was shown 
that, even though such a model can lead to early-time inflation at large values of the 
scalar curvature and late-time expansion at small values, as shown in [2Sj, it is not 
possible to have a smooth passage from the former to the latter for any model. If the 
curvature dominates the expansion there seems to be no turning back to ordinary matter 
domination, so it is reasonable to require that this should only happen at late times. 
Avoiding such an inflationary regime however, does not necessarily imply that additional 
positive powers of the curvature should not be present, but merely that they should be 
seriously constrained. In that case one can also have the usual scenario where inflation 
is driven by an inflaton field. Other constraints were derived as well, for the coefficients 
of the positive power terms, by requesting that the modified Friedmann equation does 
not differ significantly from the standard one during big bang nucleosynthesis and that 
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the correct Newtonian limit is obtained. Using the above constraints, it was shown 
that the late time evolution of the universe is not affected by the presence of additional 
positive power terms of the scalar curvature in the action. 

As a conclusion, any result derived for a model with only negative powers of the 
curvature using late time observational data, as done in m is general enough to be 
applied to a model including additional positive power terms as well. At the same 
time, however, since such tests are totally insensitive to the presence of positive power 
terms, they cannot be used to constrain them and therefore favour or disfavour their 
presence. To close I would also like to mention the following. The results of j44| seem to 
disfavour the presence of a negative power of the curvature in the gravitational action 
by showing that ACDM is the best fit model. However, the model that better fits 
the data, does not necessarily have to be the physically preferred one, unless it can 
also be theoretically motivated and justified. The ACDM model comes with a burden, 
known as the coincidence problem, and this is the motivation for creating the numerous 
alternative models present in the literature. In this sense, what we are interested in is 
not whether a model can fit the data better than the ACDM one, but if we can produce 
a model which is theoretically motivated or at least theoretically explained and at the 
same time fits the current data reasonably well. Whether this is true for f(R) theories 
of gravity in the Palatini formalism is a question that remains unanswered. 
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